Structural and Superconducting Transitions in Mgi_^Ala;B2 
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From systematic ab initio calculations for the alloy system Mgi-^, Ala;B2 , we find a strong tendency 
for the formation of a superstructure characterized by Al-rich layers. We also present a simple 
model, based on calculated energies and an estimate of the configurational entropy, which suggests 
that the alloy has two separate concentration regimes of phase separation, with critical points near 
X = 0.25 and x = 0.75. These results, together with calculations of electronic densities of states in 
several ionic arrangements, give a qualitative explanation for the observed structural instabilities, 
as well as the x-dependence of the superconducting Tc for x < 0.6. 

PACS numbers: 



The superconducting properties of MgB2 [0 remain 
the subject of intense research. Although superconduc- 
tivity in MgB2 (Tc = 39 K) appears to result from a 
phonon-mediated BCS-like interaction |^-^ the details 
of this mechanism, including the possible relevance of 
anharmonic effects multiple gaps and Fermi 

nesting 0|, are still being investigated. Studying the ef- 
fects of doping is very important, as it may not only give 
additional evidence on the origins of superconductivity in 
pure MgB2 , but is also needed to explain the observed 
structural instabilities P-]lO[] and experimental difficul- 
ties in verifying the predicted increase in Tc with Na or 
Ca substitutions 1 1 1 • 

Alloys of the form Mgi_a;Ala;B2 are the most widely 
studied experimentally of all the doped MgB2 materi- 
als |p|-p^. These systems exhibit a variety of unusual 
behavior. For example, X-ray diffraction results suggest 
that Mgi_^Al2;B2 is unstable against phase separation in 
the concentration range 0.09 < a; < 0.25 and again near 
X = 0.7 ^|-|lo|. Secondly, the superconducting transition 
shows unusual behavior as a function of x: the transition 
is broad around x — 0.25, consistent with phase separa- 
tion, then the transition temperature Tc drops sharply 
within the single phase region (0.25 < x < 0.4) §, 
but superconductivity persists (with Tc ^ 10 K) up to 
X ^ 0.7. Thirdly, a superstructure appears to form near 
X = 0.5 |P,|lo|, corresponding to Al ordering in the c di- 
rection, and possibly also at other Al concentrations. 

In this study we investigate the energetics of Al-doped 
MgB2, and their possible relation to superconductivity. 
Our particular aim is, first, to determine which struc- 
tures have the lowest energy at several concentrations, 
especially x — 0.5 and x = 0.333, and, secondly, to 
use this knowledge to shed light on the phase separa- 
tion which may occur at small and at large x, and the 
relation of these structural phase transitions to the loss 
of superconductivity with doping. While the influence 
of Al doping on superconductivity in Mgi_j;Al2;B2 has 
been discussed theoretically by several authors [pT|-p^, 



none have considered the effects of these superstructural 
transitions. 

We have carried out ab initio calculations of the total 
energy for several compositions of Mgi_a:Ala;B2 , using 
the Vienna Ab Initio Simulation Package (VASP) ||l|,|l5|, 
which employs a plane wave implementation of density 
functional theory [|6|. We used ultra-soft pseudopoten- 
tials JlTf within the generalized gradient approximation 
JlSf . For all compositions, we first arranged the ions into 
an ideal MgB2 -like structure, then relaxed the positions 
of individual ions within a computational supercell until 
the energy had converged to a chosen tolerance. At most 
X considered, we did calculations for several possible ionic 
arrangements, in an effort to determine the energetically 
favored superstructure. 

Our main numerical results are summarized in Table I. 
For X = and 1, our calculated lattice parameters, band 
structure and density of states are in very good agree- 
ment with experiment or with those calculated by other 
authors. We now discuss our results at other x, starting 
with a; = 1/3. In Fig. 0(a), we show the supercell used 
to model this composition assuming equal concentrations 
of Al ions in the different Mg layers (entry d in Table I). 
After ionic relaxation, the B ions shift from their origi- 
nal positions towards the neighboring Al atoms, as indi- 
cated by arrows. This behavior is not surprising, since 
the fully ionized Al+"^ ions carry an additional +e charge 
compared to the Mg+^ ions, thus attracting the B^ ions. 
But this ionic relaxation, since it requires altering the 
length of the strong in-plane cr-bonds formed by the B 
sp^ orbitals, is very small (^ 0.01 A) in-plane, with a 
correspondingly small energy change (~ 0.02 eV per Al). 

By contrast, if the Al atoms are assumed to completely 
fill every third Al/Mg layer [see Fig. [^(b)], the entire B 
layers shift towards the Al layers by more than 0.1 A. 
Table I shows that this relaxation reduces the energy by 
about 0.2 eV per Al atom. Hence, this layered super- 
structure is much more favorable energetically than that 
with Al ions miiformly distributed in the Mg layers. 
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Similar behavior is observed at other values of x [see 
Table I and Fig. |l|(c)]. The energies of "fully layered" 
superstructures are always lower than those of struc- 
tures in which Al is uniformly distributed in the Mg 
layers, because of this relaxation effect. The large ef- 
fect of layering is made clear in Fig. ^ where we plot 
AE[x) = Egroundix) - Eiin{x). Eground{x) IS the en- 
ergy per Mgi_xAla;B2 formula unit of the fully lay- 
ered ground state superstructure at concentration x, and 
Eiin{x) = (1 — 2;)i?MgB2 +a^£'MgAi2 is the linear interpola- 
tion. Clearly, ^E(x) is just proportional to the number 
of B layers situated between neighboring Mg and Al lay- 
ers. This behavior is reasonable, since only these B layers 
can undergo the preferential relaxation which favors the 
layered superstructure fl^ . 

Although the fully layered superstructure is always the 
lowest in energy at any a;, one can attain nearly the same 
reduction in energy by segregating the Al into various 
partially layered superstructures. We illustrate this point 
by considering x = 1/3. In the fully layered superstruc- 
ture, the Al ions occupy all the Mg sites in every third 
Mg layer [Fig. |l](b)]. In an arrangement we denote the 
"2/3" structure, the Al's fill 2/3 of the sites in every 
second layer [cf. Fig. |l|(d) and Table I]. These two struc- 
tures have nearly the same energy [GQI, which is signif- 
icantly lower than that in which the Al's are randomly 
distributed in the Mg layers. We now use this fact as 
the basis of a simple model for phase separation in these 
alloys, considering for simplicity only the regime x < 1/2. 

At a given value of x, the quantity of interest is the 
Helmholtz free energy per three-atom primitive cell of the 
MgB2 structure, which we write F{x, T) = E{x)-TS{x). 
We consider a sample that has a total of Nz layers, each 
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fully layered 
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fully layered 


3.05 
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d 
e 
f 
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no layering (eclipsed) 
no layering (staggered) 
partially layered 
fully layered 


3.05 
3.05 
3.05 
3.04 


1.105 
1.11 
1.115 
1.125 


-16.117 
-16.107 
-16.148 
-16.1512 


-0.006 
-0.001 
-0.042 
-0.067 
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fully layered (2-(-3 layers) 


3.04 


1.12 


-16.294 


-0.074 
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no layering (eclipsed) 
no layering (staggered) 
fully layered 


3.04 
3.04 
3.03 


1.095 

1.10 

1.11 


-16.459 
-16.439 
-16.512 


-0.007 
-0.0004 
-0.079 
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fully layered 


3.02 


1.105 


-16.772 


-0.068 
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pure AIB2 


3.005 


1.09 


-17.245 






layer having Ni of Mg/Al sites. We assume that x^N^ of 
these layers are Al-rich, each with Al concentration Xa, 
and (1 — Xz)Nz layers are Al-poor, with concentration 
Xm- In the regime x < 0.5, we assume for simplicity that 
Xm ~ 0, from which it follows that XaXz = x. (When we 
include Xm as a variable, in a suitably generalized free 
energy, we find that F(x, T) is minimized by Xm for 
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FIG. 1. Some of the Mgi_a;Al2,B2 superstructures studied: 
(a), (b), (c) and (d) correspond to entries d, g, k and / in 
Table I. Large circles denote Al (shaded) and Mg (open) ions; 
small open circles denote B ions, (a) Top view of a; = 1/3 
structure with no Al layering, showing both B and Al/Mg 
layers (actually separated by c/2 in the c direction). The 
in-plane displacement of the B ions by ~ 0.01 A towards 
the Al ions is indicated by arrows, (b) The ground state 
superstructure at a; = 1/3: Al ions occupy every third Mg/Al 
layer; two layers of B's are displaced towards Al's by ~ 0.1 A. 
(c) The ground state superstructure at s = 1/2. (d) Example 
of an alternative, higher-energy "2/3 structure" {x = 1/3): 
The Al ions fill 2/3 of the sites in every second Mg/Al layer. 
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TABLE I. Calculated equilibrium lattice parameters a and 
r = c/a, total energy E per Mgi_a;Ali:B2 formula unit, and 
the change in energy per formula unit Erix due to ionic relax- 
ation for several possible superstructures of Mgi_i:Ala;B2 at 
different values of x. For entries denoted "eclipsed" (as op- 
posed to "staggered") the Al ions are situated directly above 
one another in successive Al/Mg layers in the c direction. 



FIG. 2. Energy Egroundix) of the fully layered (ground 
state) structure per Mgi_j;Al2,B2 formula unit, as given in 
Table I, minus the linear interpolation Eun(x) between the 
energies of MgB2 and AIB2 , plotted versus x. The full line 
simply connects the calculated points, and the dotted line 
connects the points at 2; = 0, 0.5 and 1. 
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temperatures up to 500 K JTot.) 

We assume that the internal energy E does not depend 
on the arrangement of the Al ions within the Al-rich layer, 
but only on Xz and Xa pl[ |, and we consider only the 
region < a; < 0.5. We then make the approximation 
that 

E{Xa,Xz) = Erando-rn{x) - n{Xz){EiXa + E2XI) , (1) 

where Erari(iom(x) is the energy of random Mgi_j;Alj:B2. 
We approximate Erandomi^^) as varying linearly between 
X = and x = 0.5, as is approximately true from our 
numerical results. The second term in eq. (1) is the en- 
ergy reduction due to superstructural ordering, Eord, dis- 
cussed above. It is proportional to the fraction of B lay- 
ers n{xz) — 0.5 — \xz — 0.5| which are situated between 
Al-rich and Al-depleted layers; this effectively bounds 
the range of possible values of hy Xz < 0.5. From 
our numerical calculations, Eord is also roughly propor- 
tional to Xa- In addition, we allow for a term quadratic 
in Xa, to insure that the "fully layered structure" has 
an energy lower than that of "partially layered" struc- 
tures. This quadratic term is crucial for phase sepa- 
ration. We obtain estimates of Ei and E2 from en- 
tries d-g from the Table I: we assume that at a; = 1/3, 

^random — "^^d ~^ -^{-^z — 1/3, = 1) — Eg, 

and E{xz = 1/2, Xa = 2/3) = Ef. These relations yield 
El - 0.1 eV, E2 - 0.03 eV. 

We estimate S{x) simply as the sum of the config- 
urational entropies of the individual layers ||2^]. The 
standard expression for this entropy of one layer hav- 
ing Ni sites, of which pNg are occupied by Al ions, is 
—kBNg[p\np+ (1 —p) ln(l —p)]. Thus, the total entropy 
(per Mgi_^Alj.B2 formula unit) is estimated as 

S = -ksXz [Xa InXa + (1 - Xa) ln(l - Xa)] (2) 

We have numerically minimized the free energy 
F(xa, Xz, T) = E{xa,Xz) — TS{Xa, Xz) for fixed a; and 
T with respect to Xa, subject to XaXz — x. We call this 
resulting free energy Fmin{x,T), and define the quantity 
FHnearix,T) = (1 - 2x)i^™„ (0, T) - 2xF„„„ (0.5, T) , The 
resulting isotherms of AF = Fmin{x,T) - Fiinear{x,T) 
are plotted in Fig. ^ for several T. They have the classic 

AF, meV AF, meV 




FIG. 3. AF = Frmn(x,T) - Fu„ear{x,T), as Calculated 
from the model described in the text for (a) T = E2 ~ 350 K; 
(b) T = E2/2A5 « 140 K. The dashed line shows the common 
tangent construction which determines the composition limits 
of the two phases in the phase separated region (provided 
Erandom{x) IS linear in x). 



shape associated with phase separation: concave up for 
T > Tinst', concave down for T < Tinst- In the latter 
regime, the concentrations of the two coexisting phases 
are determined by the standard common tangent con- 
struction sketched in Figure |[ The critical temperature 
Tinst and critical concentration Xinst for phase separa- 
tion are determined as the maximum T , and correspond- 
ing X, where {d^AF/dx^)T = 0. We find Tinst = £^2/2 
{independent of Ei), and Xinst = 1/4. For our choice 
E2 = 0.03 eV, this procedure gives T^^t « 175 K. The 
experimental value of Tinst is unknown but must exceed 
the temperature at which a phase separated mixture was 
reported |p|-|l0| (presumably room temperature). But our 
estimate is obtained using an extremely simple means of 
estimating E2, and would probably be improved by a 
more elaborate calculation (moreover, our data suggest 
slight deviation of Erandom from linear behavior, favoring 
phase instability). 

For concentrations x > 0.5, a similar model could also 
be applied, probably with different parameters Ei and 
E2, leading once again to a region of phase separation 
with a critical concentration Xinst = 0.75. On the other 
hand, the upward curvature in E{x) at a; > 0.5 (cf. 
Fig. H) should oppose phase separation, decreasing the 
width of the two-phase region at a given temperature. 
This behavior once again appears to agree with experi- 
ment, as the observed two-phase region near x 0.7 is 
reported to have much smaller width at comparable T 

Finally, we discuss the observed variation of supercon- 
ducting transition temperature Tc{x) with x, based on 
these results. Our calculations confirm the suggestion 
pT|-p^ , that the decrease of Tc with increasing x is due 
primarily to a reduction in the density of states (DOS) 
near the Fermi energy. By way of illustration, we show 
in Fig. ^ our calculated Kohn-Sham DOS N{e) in pure 
MgB2 and the fully layered Mgo.5Alo.5B2 . (We have at- 
tempted to minimize structure due to spurious Van Hove 
singularities produced by the computational algorithm 
p3| , ^ by using k- meshes as fine as 35 x 35 x 35.) In- 
deed N[eF',x — 0.5) < N{eF,x = 0) as expected. The 
observed large width AT{x) of the superconducting tran- 
sition near a; = 1/4 occurs, we believe, because this con- 
centration lies in the two-phase regime. To some extent, 
DOS behaves as predicted from the rigid-band model, 
simply shifting in energy, relative to e^?, without greatly 
changing its shape. However, the slightly broader DOS 
at X — 0.5 is, we believe, a real departure from the rigid- 
band picture, and due to the increased physical unit cell 
size. 

For further insight into the occurrence of superconduc- 
tivity, we have examined the calculated band structures 
at different values of x and Xz ■ We paid special attention 
to the (T-bonding p^y bands believed to be primarily re- 
sponsible for the superconductivity. Using a rigid-band 
model for small variations in x at fixed Xz , we found that 
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for any value of Xz, these cr-bands fill at a; « 0.6. (For 
example, they are filled for structure /.) If these were 
the only occupied bands, the superconducting Tc would 
seem to vanish above x w 0.6. Since the electron-phonon 
coupling constant for those bands that remain par- 
tially filled at a: > 0.6 is very small (Ap^ ~ 0.28 |||), it 
could not produce superconductivity for T > ~ 0.01 
K, according to the Allen-Dynes formula j2^j2^. Thus, 
the persistence of a finite Tc 10 K) may result from 
some kind of coupling between electrons in the Pz and 
a bands. For example, a strong pairing interaction be- 
tween (T electrons could make it energetically favorable 
for some electrons to transfer into the Pz band. 

To summarize, the present ah initio study of 
Mgi_2;AUB2 has led to three principal findings. First, 
we find that at low temperatures a layered superstruc- 
ture is energetically preferred, not only at x = 0.5 as was 
found experimentally [^0, but also at other values of 
X. Secondly, we have described a very simple model for 
phase separation in these alloy systems. The model is 
based on a balance between the calculated ah initio en- 
ergies and configurational entropy, and leads to critical 
points at a: = 0.25 and x = 0.75, also consistent with 
experiment. Finally, we find that at x < 0.5 the exper- 
imental trends in both Tc{x) and the width Ar(a;) of 
the superconducting transition can be qualitatively in- 
terpreted in terms of the calculated x-dependent density 
of states N{eF,x), but that the interband coupling must 
be crucial in maintaining finite Tc at x > 0.6. 

We thank R. Hennig for useful discussions. This 
work was supported by NSF grant DMROl-04987, NASA 
Grant NCC8-152, and the U.-S. Israel Binational Science 
Foundation. Calculations were carried out using the fa- 
cilities of the Ohio Supercomputer Center. 
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FIG. 4. Electronic density of states A^(e, x) (per MgB2 
formula unit per spin) for a; = and for the fully layered 
superstructure at x = 0.5 [the structure shown in Fig. blc)]. 
Vertical full and dashed lines denote ef &t x = and x — 0.5. 
The two curves are lined up so that they each have a filling 
of eight valence electrons per formula unit at eF(2: = 0). 
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